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Abstract 

We show that simple parametrisations at small x of the proton structure 
function work so well in limited regions of the (x,Q 2 ) plane because they 
are approximately "self-consistent" solutions of the QCD evolution equation. 
For a class of them, we predict their Q 2 dependence and compare the result 
with experimental data. 
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1. Introduction 

Since long time non-Regge behaviour of the proton structure function F 2 (x, Q 2 ) 
in deep inelastic scattering Jl| has been discussed on the basis of the scale break 
predicted from QCD. Recent experimental data |2|, || provide further evidence for 
perturbative QCD evolution equations f|, |5| and confirm @ the scaling behaviour 
of the singlet proton structure function at small x and large Q 2 predicted in ref. @. 
This double asymptotic scaling is derived from the exact solution of the dynami- 
cal equations || in the region where nonsinglet contributions to F 2 (x,Q 2 ) can be 
neglected and the lowest moments of the splitting functions dominate. 

Empirical expressions for F 2 (x, Q 2 ), where the functional form of the perturba- 
tive double scaling is respected, will reproduce well experimental data. An example 
is the parametrisation of ref. 0, where an interesting similarity in the energy depen- 
dence between the average charged multiplicity in e + e~ collisions and the proton 
structure function has been exploited. 

Attempts to implement the unitarity condition in deep inelastic scattering || 
H [K| lead to a quite different behaviour for the proton structure function. As sug- 
gested in [|n[ double scaling admits local approximations due to the slow variation 
with x and Q 2 of the relevant variables in the asymptotic region. We agree with 
this interpretation but other reasons can be found for the agreement of these simple 
parametrisations with HERA data. 

If we consider two simple expressions for F 2 (x, Q 2 ): 

/ Tn \ c(Q 2 ) 

F 2 (x,Q 2 ) = a(Q 2 ) + b(Q 2 ){j) (1) 

and 

F 2 (x,Q 2 ) = Y,<Q 2 )(^-) 1 (2) 
and fit them at low x, for each Q 2 bin, to the experimental data @, H|, we find that 
ansatz (Q) has a lower total Xdo f 5 while the exponent c(Q 2 ) in (JJ) agrees with the 



finding of p. The good description of data obtained with the parametrisation @ 
must find an explanation in terms of the evolution equation ||. The purpose of 
this paper is to find a reason for this agreement. 

In Section 2 we will discuss new approximate solutions of the standard evolution 
equations [BJ in a form similar to the ones found in |8], [| [II]]. By limiting ourselves to 
the double logarithmic approximation, an explicit example will be given. In Section 
3 we will show that there is agreement between these solutions and an explicit fit 
to the data. Last Section is devoted to a discussion of the results. 

2. A theoretical explanation. 

Let us consider the singlet evolution equation || at leading order, when only the 
largest singularity in moment space contribute to the splitting functions. Neglecting 
non singlet contributions and defining the new variables 



«-H'?)-C- h (^) (3) 



x 



a simplified form for the problem of determining the proton structure function can 
be derived |J in the double logarithmic approximation 

FZ=^Q(Z,() (4-a) 

Q(Z, = Q{t, o) + ^ / q C d?G{t, C) (4-b) 

SO +f »f!-^,0 = o (4-0 

where Q(£, C) and are the singlet quark distribution and the gluon distri- 

bution multiplied by x, respectively. In the equations above, / is the number of 
flavours and 



with /9b = 11 - 2//3. 



As well known [ FJ| the Riemann function for the Goursat problem in eq. ( |4-c| ) is 



U{teie,0 =e 5{C - C ' ) /o[27 V / (e-0(C , -C)] (6) 

where Iq{z) is the modified Bessel function of order zero. The solution of ( [4-c| ) can 
be written in the form 

/• p i / dG\ r p i dG 
G(P) = (UG) A + j A d(U IsG +^r)+J A diU^ (7) 

where P = (£', ('), P 1 = (£ , CO, p 2 = (£', Co) and A is the point (f , Co) of the (') 
plane. 

The asymptotic form of Io(z), for large argument, leads from (^) to the double 
asymptotic scaling 0, for soft boundary conditions along AP\ and AP2. In the 
asymptotic region, starting scales xq and Qg are not important and a lower bound 
on y^C has been imposed in in order to get scaling in this variable. 



Near the origin of the (£, C) plane, eqs. ( f4-a| , f4-b| ) give a poor approximation of 



the reality, if Xq is not small enough, because the non singlet contribution becomes 
sizable. We can try to avoid in part the valence quark region by imposing on the 
experimental data the cuts: x < 0.05 and y > 0.02 as in ref. 0. Once has been 
fixed, x = 0.05, with the choice Q 2 = 6.5 GeV 2 HI § and ZEUS || data cover a 
small strip of the (£, C) plane. The maximum value of y = Q 2 /(xs) is 0.56 for HI 
and 0.9 for ZEUS, s = 90200 GeV 2 is the square of the center of mass energy of the 
lepton-proton collision. With A = 200 MeV the aforesaid bounds can be translated 
in the variables £, C 

5.1e< + e«ln(^) (8) 



and it turns out that, for three (four) flavours, the argument of the Bessel function 
in (H) reaches a maximum value of 2.8 (3.0) for ZEUS and 2.6 (2.8) for HI. 

In all the above region the series expansion for Io(z), truncated to the first three 
terms, provides an approximation better than the standard asymptotic expansion 
for all z, up to z ~ 2.8. However, already at z — 4, the relative errors become 20 
and 3.6 percent, respectively. This can explain both the limitations of this trun- 
cation and the interest for a local parametrisation whose success is doomed to an 
abrupt failure. Outside the region outlined above, deviations from QCD perturba- 
tive evolution are in fact under control in the simple expression for F 2 (x,Q 2 ) we 
will obtain. The onset of new dynamical effects, like damping by screening cor- 



rections [13], can be easily detected and, finally, a reasonable and simple input for 



evolution is suggested, just at the border of the existing published data. 

We return now to eq.(|7|). Boundary conditions can be given along the lines 
Co = and £ « 2.75 (3.0) for HI (ZEUS). This choice satisfies the cut (H) and is 
equivalent to a shift in the starting point x$ as far as the evolution is concerned. A 
more refined solution, where part of the boundary is given by eq.(|8]) with y = 0.02, 
would meet with difficulties in the choice of the initial distributions. Since the series 
for Io(z) will be truncated, Iq(2*/z) rs 1 + z + z 2 /4, the evaluation of for 
given boundary conditions, is immediate. If the starting conditions are 

G(e,0)=a 1 + A(e-^ ) (9) 

and 

G(Z ,0=a l + a 2 ( (10) 

with a n-terms approximation for Iq{z), we end up with a polynomial in £, of 
degree (n+1), for G(£, (). We argued that a three terms expansion for Iq(z) is 



needed to reproduce correctly the Q 2 evolution, then G(£, (), and hence F 2 , should 
be a polynomial of fourth degree in the variable £ with coefficients depending on (. 

In order to reduce the number of free parameters we can perform the last integral 
in eq.(0) obtaining 



a— J^T^iiiWcce'-eo) 



7 V C 

and then approximate the integral in eq. ( |4-b|) in the form 

K' e -*C 



'o vC 

7 ye^ g(i - + - ^'K' + i))) (ID 

Errors in eq.flTTD have been checked numerically and are compatible with the ap- 
proximation adopted. At high Q 2 , the result will lose in accuracy but, with this 
expedient, we can utilize in the fit the Q 2 bins where few data appear. 

If we rewrite eq. (fj) in the form 



^ = ^I>(C)(£-£o) 



;i2) 



i=0 



we get the relations (c« = Cj(0)) 



7 4 / 

Cl(C) = Cl+ 9^ 



i 2 f c 

Co(C) = c + — («iC + Oi 2 —) 



(13-a) 
(13-b) 



c 2 (C) = c 2 + 



TV 
18<5 3 



3a 2 5 2 A S 



-2ai + — + 



S i 1 



-<5C 



+ 



(13-c) 
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Boundary conditions are imposed by giving c (C) (hence c , cti and a 2 are 
known) and the constants c« (i = 1,2). Finally the constraint 

a=i(0 



9 

C=o J 



9C 

determines A if the derivative of ci(£), for £ = 0, is known from experiment. The 
application of the above condition is difficult since experimental data have large 



errors. However gluon density must satisfy other approximate constraints [ I I 



15| , jig] and, since it is universal, its behaviour can be inferred also from vector 
meson [T7| and diffractive jet [T^j production. 



We notice that the initial conditions @ and (|10D must be considered as an 
example since (|9|) should be multiplied by a factor (1 — x) a , that is (1 — x e~^) a , 
and ( |10D is perhaps too naive. As a consequence of the neglect of the term (1 — x) a , 
predictions for F2(x,Q 2 ) will be above the experimental data for the points with 
large x. The effect of this approximation will be in fact more important than the 
neglect of the valence quarks contribution in the selected region. 

3. Comparison with a fit to the data. 

We argued that eqs.([L^) and ( |1 3- a| - |13^ ) will reproduce the experimental proton 



structure function in the region of the (x, Q 2 ) plane specified above. To substantiate 
this belief we performed a fit of the experimental data, for each Q 2 bin separately, 
to the functional form (Il2|). 



To facilitate the comparison with the theoretical predictions, £ m (12) has been 
chosen as £o = ln(0. 05/0. 0032) for HI data that will be considered first. The 
total Xd.o.f. f° r au the 16 values of Q 2 , from 6.5 GeV 2 to 500 GeV 2 , is 3.9 using 
the full errors, with a mean value of 0.24 for each Q 2 bin fitted. Parametrisation 
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(|T|) would give instead a total Xd.o.f. °f 4-7 This finding justifies our claim in the 
Introduction that parametrisation (||) is preferable to ([!]). 

The result from the fit is displayed in fig.l, where the values of the coefficients 
in eq.fllj) 

4(0 = ^,(0 (< = 0,1,2) (14) 

have been plotted with closed circles as function of Q 2 . The continuous curves in 
fig.l are obtained from eq. (|13-a - 13^) with the following boundary conditions 



c = 2.24, ci = 0.46, c 2 = 0.063 

a x = 5.17, a 2 = 9.83, A = 2.53 (15) 
and / = 4. While the parameters a.\ and ot 2 have been fixed from a fit to F 2 (x = 



0.0032, Q 2 ), eq.( |13-a|) , a guess for A has been obtained from the experimental gluon 



momentum densities || [16 



In fig. 2 the continuous curves show the gluon momentum density as function of 
x, obtained with the parameters in (|T5|) , at Q 2 = 6.5 GeV 2 and evolved at Q 2 = 
20GeV 2 , without any approximation on the last integral in eq.(^|). Comparison 
with the proposal in P, |16| shows that the slope A is almost correct, while the 
value of «i is below what one would expect by a 20 percent. This discrepancy 
can be ascribed to the difficulty of fitting the proton structure function, at fixed 
x, as function of Q 2 . Both the simplified input and the large errors in the data 
contribute to the difference that can be quantified in a leading-order calculation. 
The dashed curve in fig. 2 represents the result for G(x,Q 2 ) obtained from eq.(|12D 



with the leading order method of ref . |14| . Next-to- leading order calculation would 



lower this prediction, but the value for aj., suggested by the fit, seems really too 



small. 

The same fit has been repeated for the ZEUS data || . The total Xdo f increases, 
with a mean value of 0.8 for each Q 2 bin fitted. A (() and Ai(() appear as in 
fig.l, but the fit for A 2 (() deteriorates somewhat. While the spread of the fitted 
points increases, the curves defined in eqs.( |13-a| - |13-cD represent always well the 



general behaviour. The different regions of the (£, () plane covered by the two 
sets of data account only partly for the difference in the two fits. However, other 
parametrisations presented us with the same problem. 



2 \ 



Going back to the HI data, examples of the proton structure function ^(x, Q 
obtained with this method, are given in fig. 3 with continuous curves, at different 
Q 2 values, as function of the variable ^fs7 q introduced in ref. [0], 

2Qi 



X 



The choice of this variable, where Q\ = 270 MeV, allows for a comparison of the 
energy dependence of the mean charged multiplicity with the x behaviour of F 2 . 



Quadratic fits in In s to the multiplicity data for p— p scattering are rather old [IS 



20] and their derivation in the framework of a Regge model has been suggested [21 



Till 1992 the phenomeno logical fit of ref. [20| has been shown in the Review of 



Particle Properties [[22| and, from then, an analogous version adapted to e + — e 



annihilation became available . 

In fig. 3 we have superimposed to the data also the predictions for the mean 
charged multiplicity of ref. |20| (dashed line) and of ref. |23| (dotted line) nor- 
malized to the F% data at each Q 2 value. The idea of similarity between the two 
observables, advanced in ref. J7|, seems confirmed in this more qualitative approach. 
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Remarkably, this similarity could be process independent at high energy. Notice 
that [^| the empirical logarithmic solution has a better Xdo f than the leading-log 
approximation proposed in |24j] and provides a fit to the mean charged multiplic- 



ity comparable to the modified leading-log approximation [g5j in the energy range 
under consideration. 

4. Conclusions 

Motivations for the success of simple parametrisations for the proton structure 
function have been found in the structure of the QCD evolution equations. In this 
paper we have shown that a region of the (x, Q 2 ) plane exists where both asymptotic 
and truncated series expansions give the answer with a comparable accuracy. In the 
double logarithmic approximation, we are considering, this region encloses almost 
all HERA data while, at smaller x, the asymptotic solution would certainly be the 
only appropriate. 

Except for an exact next-to-leading order computer calculation, the logarith- 
mic fit appears to be at the level of other simple parametrisations, at least in the 
explored domain of the relevant variables. This phenomenological finding has a cor- 
respondence in the theoretical picture of deep inelastic scattering, where F 2 (x,Q 2 ) 
is obtained from the discontinuity of the Pomeron exchanged between the proton 
and a quark loop coupled to the photon. Representing the Pomeron as a gluonic 
ladder, the Born term, correponding to the exchange of two gluons, would produce 
roughly an x independent asymptotic structure function. Every rung added to the 



ladder brings about a term increasing like ln(l/x) when x tends to zero | 26fl . Hence, 
a two-rung gluonic ladder is suggested from our logarithmic fit. 

Two important points must however be noticed. First, at the border of the 



(£, C) region studied, this approach shows his limitations and the fit deteriorates 
somewhat. The second point regards the possibility to extend the range of validity 
of this method by keeping more terms in the expansion for the Bessel function. The 



answer is negative. Mathematical tables [27| show that, in the double logarithmic 
approximation, double scaling becomes asymptotically ineluctable. 
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Figure Captions 

fig. 1: Coefficients of the logarithmic fit, plotted as function of Q 2 for the data 
in [0], the continuous line shows the prediction of the model (see text). The 
point at 500 GeV 2 for Ax has been shifted to the right for the sake of clarity. 

fig. 2: The gluon momentum density in the model, as function of x, at Q 2 = 
6.5 GeV 2 and evolved at Q 2 = 20 GeV 2 (continuous line). Dashed line shows 
the result of a leading order calculation [[Hfl . 

fig. 3: Proton structure function as function of y/s^ (see text) at different Q 2 
values, data are from |2j and continuous lines from our model. Predictions 
obtained from the mean charged multiplicity, normalized to F2(x,Q 2 ), are 
from ]50] (dashed line) and from (dotted line). 
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